Mon. Not. R. Astron. Soc. OOP. [HM1 ('2005) Printed 5 February 2008 (MN 1^1^ style file v2.2) 

A relativistic model of the radio jets in NGC 315 



J.R. Canvin , R.A. Laing^, A.H.Bridle^, W.D. Cotton^ 

School of Physics, University of Sydney, A28, Sydney, NSW 2006, Australia 
^ University of Oxford, Department of Astrophysics, Denys Wilkinson Building, Kehle Road, Oxford OXl 3RH 
^ European Southern Observatory, Karl-Schwarzschild-Strafle 2, D-85748 Garching-bei-Miinchen, Germany 
* National Radio Astronomy Observatory, 520 Edgemont Road, Charlottesville, VA 22903-2475, U.S.A. 



Received 



ABSTRACT 

We apply our intrinsically symmetrical, decelerating relativistic jet model to deep VLA 
imaging of the inner ±70 arcsec of the giant low-luminosity radio galaxy NGC 315. An 
optimized model accurately fits the data in both total intensity and linear polarization. 
We infer that the velocity, emissivity and field structure in NGC 315 are very similar 
to those of the other low-luminosity sources we have modelled, but that all of the 
physical scales are larger by a factor of about 5. We derive an inclination to the line 
of sight of 38° ± 2° for the jets. Where they first brighten, their on-axis velocity is 
(3 = V / c ~ 0.9. They decelerate to (3 ~ 0.4 between 8 and 18kpc from the nucleus 
and the velocity thereafter remains constant. The speed at the edge of the jet is «0.6 
of the on-axis value where it is best constrained, but the transverse velocity profile 
may deviate systematically from the Gaussian form we assume. The proper emissivity 
profile is split into three power-law regions separated by shorter transition zones. In 
the first of these, at «3 kpc (the fiaring point) the jets expand rapidly at constant 
emissivity, leading to a large increase in the observed brightness on the approaching 
side. At «10 kpc, the emissivity drops abruptly by a factor of 2. Where the jets are well 
resolved their rest-frame emission is centre-brightened. The magnetic field is modelled 
as random on small scales but anisotropic and we rule out a globally ordered helical 
configuration. To a first approximation, the field evolves from a mixture of longitudinal 
and toroidal components to predominantly toroidal, but it also shows variations in 
structure along and across the jets, with a significant radial component in places. 
Simple adiabatic models fail to fit the emissivity variations. 
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- MHD 



1 INTRODUCTION 

It was first recognised by iFanarofI fc RilevI il974l) that 
the ways in which extragalactic jets dissipate energy to 
produce observable radiation differ for high- and low- 
luminosity sources. Jets in low-luminosity (FRI) sources 
are bright close to the nucleus of the parent galaxy, 
whereas those in powerful (FRII) sources are relatively 
faint until their terminal hot-spots. It rapidly became 
accepted that FRI jets m ust decelerat e by entrainment 
of t he surrounding IGM (Baan 1980: BcEcIman "1982'; 
i Bicknell 1984. 1986.:. De Young 1996. 2004; R osen ot al. 1999; 
iRosen fc Hardeell2000D o r by injection o f mass lost by stars 
within the jet volume ([Phimievl Il983t iKomissarovl Il994l: 
iBowman. Leahy fc Komiss£TO^^96h . 

More recently, evidence has accumulated that FR I jets 
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are initially relativistic and decelerate on kpc scales. FRI 
sources are thought to be the side-on counterparts of BL 
Lac objects, in w hich relativistic motion on parsec scales 
is well-established JUrrv fc Pa dovani" 1995V Proper motions 
corresponding to speeds comparable with and in some cases 
exceeding c have b een measured on millia rcsecond scales 
in several FRI jets iGiovann in^^L||20^|) and on arcsec- 
ond scales in M 87 and Cen A jBireteil Zhou fc Ow'enlll995l: 
iHardcastle et alJl200j) . In FRI sources (as in FRIIs), the 
lobe containing the main (brighter) jet is less d epolarized 
than the counter-jet lobe iMorganti et al.lll997^ . This can 
be explained as an effect of Faraday rotation in the sur- 
rounding halo of hot plasma if the main jet points toward 
the observer, suggesting that the brightn ess asymmetry is 
caused by Doppler beaming jLaing|ll988l) . The decrease of 
this a symmetry with distance from the nucleus iLaing et alJ 
I1999I) indicates deceleration. 

The present paper is the third in a series devoted to 
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modelling the jets in FRI radio galaxies. We assume that 
the jets are intrinsically symmetrical, axisymmetric, rela- 
tivistic flows and we parameterize their geometry and the 
three-dimensional variations of their velocity, emissivity and 
magnetic-field structure. We then compute the brightness 
distributions in total intensity and linear polarization. By 
fitting simultaneously to deep, high-resolution radio images 
in Stokes I, Q and U, we can optimize the model parameters. 
The fits are empirical, and make as few assumptions as pos- 
sible about the (poorly-known) internal physics of the jets. 
The technique was originally developed to model the radio 
galaxy 3C31 bylLaing & Bridle (2002a, hereafter LB) and 
was then slight ly revised and applied to B2 0326-1-39 and 
B2 1553-^24 bv ICanvin fc Laind i2004 hereafter CL). We 
now present a model for the jets in the well-known nearby 
radio galaxy NG C 315. Th is source has a huge angular and 
physical size (Bri dle et alj ll976) and, as we shall show, the 
angular scales on which fiaring, recollimation and decelera- 
tion of the jets take place are much larger than in the sources 
we have studied so far. Our results for NGC 315 therefore 
give a more detailed picture of the initial propagation of an 
FRI jet. 

In order to improve on our empirical models, we need to 
understand the energy gain and loss processes affecting the 
ultrarelativistic particles which produce the observed syn- 
chrotron radiation. A self-consistent, axisymmetric adiabatic 
model fails to fit either the total intensity or the polarization 
distributions in the je ts of 3C 31 within 5 kpc of the nucleus 
(iLaing fc Bridlell2004) . This suggests that injection of rela- 
tivistic particles and/or amplification of the magnetic field 
are required, which is not surprising in vie w of the detection 
of co spatial X-ray synchrotron emission iHardcastle et af] 
l2002ll . Further out in 3C31, the adiabatic model gives a 
tolerable fit. A simple, quasi-one-dimensional analysis (LB, 
CL) is adequate to assess whether fitting of more elaborate 
models is worthwhile, and we apply this to NGC 315. 

Given a kinematic model for the jets and estimates of 
the external density and pressure from X-ray observations, 
we can apply conservation of mass, momentum and energy 
to deduce the variations of internal pressure, density, Mach 
number and entrainmen t rate with distance from the nucleus 
(iLaing fc Bridlell2002b^ . 

In Section |5| we introduce NGC 315 and briefly sum- 
marize the VLA observations. Our modelling technique is 
outlined in Section |3 emphasizing the (small) differences 
from the earlier work of CL. The observed and model bright- 
ness distributions are compared in Section^Jand the derived 
geometry, velocity, emissivity and fleld distributions are pre- 
sented in Section|K| We summarize our conclusions and out- 
line our future programme in Section 15] 

We adopt a concordance cosmology with Hubble con- 
stant, Ho = 70kms"^Mpc"\ Qa = 0.7 and Q,m = 0.3, 
although only the choice of Hq is significant at the distance 
of NGC 315. 

2 OBSERVATIONS AND IMAGES 
2.1 NGC 315 

Our modelling technique requires that both radio jets are: 
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Figure 1. A grey-scale of the large-scale structure of the jets in 
NGC 315 at a resolution of 5.5arcsec FWHM. The image is the 
mean of four at different frequencies between 1.365 and 1.665 GHz 
and is described in Laing et al. (in preparation). The grey-scale 
range is - 2.5 mjy (beam area)~^ as shown by the labelled wedge 
and the modelled area is indicated by the box. 

(i) detectable with good signal-to-noise ratio in total in- 
tensity and linear polarization; 

(ii) straight and antiparallel; 

(iii) separable from any surrounding lobe emission and 

(iv) asymmetric, in the sense that their jet/counter-jet 
intensity ratio is significantly larger than unity over a sig- 
nificant area. 

The nearby giant elliptical galaxy NG C 315, whose large- 
scale radio structure was first imaged bv lBridle et al.l |I197EI) , 
is one of the brightest sources satisfying these criteria. The 
galaxy has a redshift of 0.01648 ( Tragor et al. 2000), giv- 
ing a scale of 0.335 kpc arcsec"^ for our adopted cosmol- 
ogy. The overall extent of the radio source is approximately 
3500 arcsec, or about 1200 kpc in projection, but the area to 
be modelled (see Fig. is limited in extent by the slight 
bends in the jet at roughly 70 arcsec (23 kpc in projection) 
from the nucleus (see Section f3.2^ . 

High-resolutio n images of t he jets on kpc scales were 
presented by Bridle et alJ (ll979h . lFomalont et alJ(ll980h an d 
'Vcnturi ct al. ( 199 j) and on pc sca l es bv iLinfie ld' l'l98l[l 
(Venturi et_aL (1993), ICotton et all (ll999^ and X u et all 
l20onir More comprehensive references to radio, optical and 
X-ray observations of NGC 315 are given by Laing et al. (in 
preparation). 

2.2 Observations and data reduction 

The observations and their reduction are also described in 
detail by Laing et al. (in preparation). We fit to images 
made using the 5 GHz dataset from that paper at 2.35 and 
0.40 arcsec FWHM [the higher-resolution images are also 
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Table 1. Image resolutions and noise levels, aj is the off-source 
noise level on the I image; ap the average of the noise levels for 
Q and U . 



FWHM 


rms noise level 


(arcsec) 


[/xjy (beam area)"-'] 




CTj crp 


2.35 


10.0 7.5 


0.40 


12.5 6.8 



discussed and compared with Chandra observations by Wor- 
rall et al. (in preparation)] . The dataset includes long obser- 
vations in all four configurations of the VLA and provides ex- 
cellent sampling of spatial scales between 0.4 and 300 arcsec. 
The resolutions and noise levels are given in Table Q 



The accuracy of linear polarization measurements is 
critical to our technique. The E-vector position angles have 
been corrected to zero wavelength using an image of rotation 
measure derived from 5-frequency observations at a resolu- 
tion of 5.5 arcsec FWHM between 1.365 and 5 GHz (Laing 
et al., in preparation). Significant corrections are required 
for the mean Faraday rotation (— 76radm~'^) and a linear 
gradient along the jet, both probably of Galactic origin. The 
fitting errors are typically <2.5radm~^ in the area of inter- 
est, corresponding to position angle rotation <0.5° at 5 GHz. 
The observed fluctuations of rotation measure on scales of a 
few beamwidths are very small (~1 - 2 rad m~^) and residual 
depolarization is undetectable, so errors in E-vector position 
angle and degree of polarization p — {Q^ + U^Y^'^/I due 
to changes in Faraday rotation across the 5.5 arcsec beam 
should be negligible. Observations at two centre frequen- 
cies (4.860 and 4.985 GHz) were combined to make the fi- 
nal dataset (Laing et al., in preparation) but the error in 
E-vector position angle due to this efi'ect is <1°. Depolar- 
ization due to rotation across the maximum bandwidth of 
50 MHz is also negligible. T he modelling described lat er is 
not affected by Ricean bias dWardle fc KronberdllQT^) . as 
we fit to Stokes /, Q and U directly, but our plots of the de- 
gree of polarization include a first-order correction for this 
effect. We plot the direction of the apparent magnetic Beld, 
rotated from the zero-wavelength E-vector position angles 
by 90°. 



3 THE MODEL 
3.1 Assumptions 

Our fundamental assumptions are those of LB and CL: 

(i) The jets may be modelled as intrinsically symmet- 
rical, antiparallel, axisymmetric, stationary, laminar fiows. 
The real flow is clearly much more complicated, but all that 
is necessary for our technique to work is that an average 
over a sufficiently large volume results in statistically iden- 
tical rest-frame emission from the main and counter-jets. In 
poorly resolved cases, some care is required to distinguish 
between intrinsic variations of physical parameters and local 
fiuctuations such as knots and filaments, as we will discuss. 



(ii) The jets contain relativistic particles with an energy 
spectrum n{E)dE = noE~''^"'^^'^ dE (corresponding to a fre- 
quency spectral index a) with an isotropic pitch-angle distri- 
bution. The maximum degree of linear polarization is then 
po — (3q -f 3)/(3a -|- 5). We use a — 0.55, the mean spec- 
tral index for the modelled region between 1.365 and 5 GHz 
(Laing et al., in preparation). 

(iii) The magnetic field is tangled on small scales, but 
anisotropic (the reas ons for taking the field to be of this 
form are discussed bv lLain j|l98ll . lB'egelman et alJll984l and 
LB; see also Section f5.5|l . 

(iv) The emission is optically thin everywhere. The core, 
which is unresolved in our VLA observations and partially 
optically thick, is not modelled. It is represented as a point 
source with the correct flux density. 

We define /3 = v/c, where v is the flow velocity. F = (1 — 
^2^-1/2 ^YiQ bulk Lorentz factor and is the angle between 
the jet axis and the line of sight. 



3.2 Outline of tlie method 

The jet/counter-jet intensity ratio for intrinsically symmet- 
rical, relativistic jets depends on the product /3 cos 9, so the 
angle to the line of sight and velocity cannot be decoupled 
by observations of / alone. The basis of our technique is that 
relativistic aberration causes the approaching and receding 
jets to appear different not only in total intensity, but also 
in linear polarization. We use these polarization differences 
as independent constraints in order to break the degeneracy 
between /3 and 9. An outline of our procedure is as follows. 

(i) Develop a parameterized description of the geometry, 
velocity, emissivity and magnetic-field ordering. 

(ii) Calculate the synchrotron emission in /, Q and U 
from the model jets by numerical integration, taking account 
of relativistic aberration. 

(iii) Convolve and compare with deep VLA images, using 
to measure the quality of the fit. 

(iv) Optimize the model parameters using the downhill 
simplex method. 

The details of the calculations are fiiUy described by LB and 
CL. 

The assumption of axisymmetry is clearly violated in 
NGC315 at «70 arcsec from the nucleus. Both jets bend 
clockwise by «5° (in projection) at this distance, but we 
do not know the magnitude of any associated bends in the 
orthogonal direction. Our modelling could be extended to 
larger distances provided that the jets are indeed antipar- 
allel and intrinsically identical after the bends. This seems 
plausible from their appearance in projection and we could 
fit independently for the angle to the line of sight after the 
bend if necessary. The brightness and polarization structures 
of the jets are qualitatively consistent with an extrapolation 
of the fitted model described below, so it is likely that this 
approach would succeed. It adds significant complexity to 
the modelling procedure, however, and also requires careful 
justification of the assumption of intrinsic symmetry at large 
distances. For these reasons, we defer this analysis to a later 
paper. 
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Figure 2. A sketch of the assumed jet geometry, defining the 
quantities goi ''0 and x'o used in Table ISl For NGC315, the mod- 
elled area (indicated schematically by shading) is entirely within 
the flaring region. The full line indicates the edge of the jet (s = f ) 
and the dashed line the s = 0.5 streamline. 

3.3 Functional forms for geometry, velocity, 
magnetic field and emissivity 

The functional forms used to parameterize the geometry, ve- 
locity, emissivity and magnetic field are very similar to those 
described by CL, which are in turn an evolution from those 
introduced by LB. The earlier papers discuss the motivation 
of these expressions in detail; here we just summarize their 
forms for completeness and to highlight a few differences 
from our earlier work. 

3. 3. 1 Geometry 

The assumed geometry, sketched in Fig.|5|is derived from a 
fit to the outer isophote of the jet emission. It is identical to 
that used by CL, but in this case we model only the flaring 
region where the jet expands rapidly and then recoUimates. 
NGC 315 also shows a conical outer region of constant (very 
small) opening angle, as seen in other FRI jets, but this is 
beyond the bends which limit our ability to fit an axisym- 
metric model directly to the observed data. 

We define z to be the distance from the nucleus along 
the axis and to be the half-opening angle of the jet in the 
outer region. As in CL, our assumed flow streamlines are 
a family of curves one of whose members is deflned by the 
outer boundary of the jet. We use a streamline coordinate 
system (p, s) where the streamline index s is constant for a 
given streamline and p increases monotonically with distance 
along it. In the outer region, the streamlines are assumed to 
be radial from a point on the axis with z = —A, where 
A = xq/ sin 5o — '"0- For a streamline which makes an angle ^ 
with the jet axis in the outer region, we define s = ^/^o, so 
s = on-axis and s = 1 at the edge (Fig. |2l • The distance 
of a streamline from the jet axis in the fiaring region is: 

x{z,s) = 02(3)2^-1-03(8)2^ 



and varies monotonically along a streamline from at the 
nucleus to ro at the spherical surface defining the end of the 
fiaring region. 

3.3.2 Velocity 

The form of the velocity field /3(p, s) is the same as that 
used by CL. The on-axis velocity profile is divided into three 
parts: 

(i) roughly constant, with a high velocity /3i, close to the 
nucleus; 

(ii) a linear decrease and 

(iii) roughly constant but with a low velocity /3o at large 
distances. 

[the velocities in (i) and (iii) are not quite constant because 
transitions between regions are smoothed to avoid numeri- 
cal problems]. The profile is defined by four free parameters: 
the distances of the two boundaries separating the three re- 
gions, pvi and pvQ, and the characteristic inner and outer 
velocities /?i and /3o- Off-axis, the velocity is calculated us- 
ing the same expressions but with inner and outer velocities 
/3i exp(— ln«i) and /So exp(— In no), i.e. with a truncated 
Gaussian transverse profile falling to fractional velocities vi 
and vo at the edge of the jet in the inner and outer regions, 
respectively. 

The full functional forms for the velocity field (3{p,s) 
are given in Table |5| 



3.3.3 Magnetic field 

We define the rms components of the magnetic field 
to be (Bf)^^^ (longitudinal, parallel to a streamline), 
{-Br )^^^ (radial, orthogonal to the streamline and outwards 
from the jet axis) and (B^)^^^ (toroidal, orthogonal to 
the streamline in an azimuthal direction). The magnetic- 
field structure is parameterized by the ratio of rms ra- 
dial/toroidal field, j(P.s) = {B^y^'^/{B^y^^ and the lon- 
gitudinal/toroidal ratio k{p,s) = {BfY^'^ / {B^Y^'^ . The pa- 
rameterization is similar to that used by CL, but we also 
allow variation of the field component ratios with s (cf. LB) 
in order to improve the fit to the variation of polarization 
across the jets, which is better resolved in NGC 315 than in 
the sources modelled by CL. For each of the axial (s = 0) 
and edge (s = 1) streamlines, the field ratios have constant 
values for p < pB^ and p > pbq , with linear interpolation be- 
tween them for pBi ^ p ^ PBq • For intermediate streamlines, 
we then interpolate linearly between the ratios for s — and 
s = 1. 

The functional forms assumed for the field ratios are 
again given in Table |5| 



The curve is constrained to fit the jet boundary for s — 1 
and the coefficients are determined by the conditions that 
x{s) and its first derivative are continuous at the boundary 
between the regions, p is defined by: 

_ 2ro 

^ " (ro + A) cos(5os) - A 



3.3.4 Emissivity 

We write the proper emissivity as e{p, s)h{p, s), where e is 
the emissivity in / for a magnetic field B = {Bf + B^ + 
Bt)^^^ perpendicular to the line of sight, h depends on field 
geometry: for 7, ^ /i ^ 1 and for Q and U —po ^ h ^ +Po. 
We refer to e, loosely, as 'the emissivity'. For a given spectral 
index, it is a function only of the rms total magnetic field and 
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Figure 3. A comparison of the model and data in total intensity at 2.35 arcsec resolution for the inner 66.5 arcsec of each jet. (a) model 
contours; (b) observed contours. The contour levels in panels (a) and (b) are -2, -1, 1, 2, 4, 8, 12, 20, 36, 60, 100, 160, 250, 400, 600, 
1000 X 30/iJy(beam area)^^. (c) and (d) / profiles along the jet axis for the data (solid line) and model (dashed line), (e) and (f) 
Grey-scales of sidedness ratio, obtained by dividing total-intensity images by copies of themselves rotated by 180° about the core. The 
ratio is in the sense main jet / counter-jet. (e) model; (f) data, (g) profile of sidedness ratio along the jet axis for the data (solid line) 
and model (dashed line), (h) - (j) as (e) - (g), but with smaller grey-scale and profile ranges to emphasise the variation of sidedness in 
the outer parts of the modelled region. 
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the normalizing constant of the particle energy distribution, 
eocnoB^+". 

The on-axis emissivity profile consists of five regions, 
each with a power-law variation of e with p. The profile 
for NGC 315 is consistent with a continuous variation of 
emissivity, so we enforce continuity everywhere (<? = 1 in 
the notation of CL). We use the expressions given by CL, 
but change the notation slightly in the interests of greater 
clarity. The profile is defined by five power-law indices and 
four boundary positions. As explained in more detail in Sec- 
tion |^|1J there are three primary emissivity regions with in- 
dices _Ein, -Eniid and iJout- The remaining two regions are 
introduced to model the transitions between them. The first 
(slope -Eriso) replaces the discontinuous increase in emissivity 
used to model the initial brightening of the jets in 3C 31 and 
B2 0326-1-39 (LB, CL) because the equivalent structure in 
NGC 315 is better resolved. The second (slope iSfaii), which 
models an abrupt decrease in emissivity, is equivalent to £3 
in the models of B2 0326-^39 and B2 1553-^24 (CL). 

Off-axis, the profile is multiplied by a factor 
exp[— In e(p)], so that e(p) is the fractional value of the 
emissivity at the jet edge. e(p) has a constant value eo for 
p > P03 and varies linearly through regions 2 and 3 from ei 
at Pd to eo at pe^ . For p < pei , the jet is too narrow for our 
data to constrain any transverse profile and we set e(p) — 1. 

The full description of the emissivity distribution e(p, s) 
is given in Tabled 

3.4 Fitting, model parameters and errors 

The noise level for the calculation of for I is estimated 
as l/\/2 times the rms of the difference between the im- 
age and a copy of itself reflected across the jet axis. For 
linear polarization, the same level is used for both Q and 
U. This is the mean of l/\/2 times the rms of the dif- 
ference image for Q and the summed image for U, since 
the latter is antisymmetric under reflection for an axisym- 
metric model. This prescription is identical to that used by 
LB and CL. The region immediately around the core is ex- 
cluded from the fit. is calculated using the high-resolution 
(0.4arcsec) images where they have adequate signal-to-noise 
ratio (0.9 - 17arcsec from the nucleus) and the lower- 
resolution (2.35arcsec) images elsewhere (17 - 67arcsec). 
The same fitting regions are used for total intensity and lin- 
ear polarization and x^ values for /, Q and U are summed 
over grids chosen so that all points are independent. 

We derive rough uncertainties, as in LB and CL, by 
varying individual parameters until the increase in cor- 
responds to the formal 99% confidence level for indepen- 
dent Gaussian errors. These estimates are crude (they ne- 
glect coupling between parameters), but in practice give a 
good representation of the range of qualitatively reasonable 
models. The number of independent points (1414 in each of 3 
Stokes parameters) is sufficiently large that we are confident 
in the main features of the model. 

4 COMPARISON BETWEEN MODELS AND 
DATA 

The quality of the fits is extremely good and the reduced 
X^ ~ 1.27 for 4242 independent points. 



(a) Total intensity: model 




5 10 15 20 

Distance from nucleus (arcsec) 

Figure 4. A comparison of the model and data in total intensity 
at 0.4 arcsec resolution for the inner 21 arcsec of the main jet. 
(a) model contours; (b) observed contours. The contour levels in 
panels (a) and (b) are -2, -1, 1, 2, 4, 8, 12, 20, 36, 60, 100, 160, 
250, 400, 600, 1000 x 30 pjy/beam area, (c) / profile along the 
jet axis for the data (solid line) and model (dashed line). 



4.1 Total intensity 

The observed and modelled total intensities and jet/counter- 
jet sidedness ratios are displayed in Figs|3]-|7] Fig.|3]shows 
the inner 66.5 arcsec of both jets at a resolution of 2.35 arcsec 
FWHM and Fig.^shows the inner 21 arcsec of the main jet 
at 0.4 arcsec FWHM. Fig.|S]shows grey-scales of the observed 
total intensity for the inner 30 arcsec of both jets at the same 
resolution, but with levels chosen to emphasize the fine-scale 
structure. These are compared with contours of the model 
I. Figs |S| and |7| show the average transverse profiles of total 
intensity and jet/counter-jet sidedness ratio over areas where 
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Table 2. Functional forms of the velocity /3, emissivity e, radial/toroidal and longitudinal/toroidal magnetic-field ratios j and k in the 
streamline coordinate system (p, s). Column 4 lists the parameters which may be optimized, for comparison with Table El 



Quantity Functional form 



Range Free parameters 



/3icxp( — s Inui ) — /3o cxp( — s Ihuq) 
10 



C2 + C3P 



/3i cxp( — s In vi) — f!iQ cyip( — s In^Jo) 
10 



Velocity field" 

exp[ci (p — pvi )] p < Pvi Distances pvi , Pvq 

Pvi ^ P ^ Pvq Velocities /3i, /3o 

exp [c4 (pvQ — p)] P > Pvq ^ Fractional edge velocities i)i , t)o 



Emissivity 



(p<=2 ) 



exp 



exp 



2i.(e, + (eo-eo(^)) 



-Efall 

"^2 ( ^ exp(-s2 Ineo) 



da 



p < pci Distances Pd , Pcj , P03 , pc. 



Pd < p ^ Pc2 Indices Bin, Ehsc, £^mid, £^fall. -^ou 



P02 < P P03 Fractional edge emissivities ei , eo 



Pea < P ^ Pe4 



P > Pc, 



Radial/toroidal field ratio 



jip,s) j™= + g(-dgc_^.axiB) 



where 



^■axis ^'axis i /^'axis ^'axisN / ^ PB]^ \ 

J -Ji +Uo Ji ^ Vpbo-pbi ^ 

•axis-) P-PBi \ 
•'1 ^^^PBQ-PBiy 



P < PBi Distances pB^ , PBq 

PBi < P < PBo Ratios if ^.cdgc_ 



P ^ PBo 



k{p, s) 



^jaxis _,_ ^(fccdgo _ ^axis) 
where fc^-xis _ ^axis _j_ (^axis _ 
fcodgo ^ fcOdgc ^ ^^odgc 



Longitudinal/toroidal field ratio 



1 ' \PB„-PBI J 
^cdgC) / P-PB, \ 
1 ^\^PBo-PBiy 



p ^ PBi Ratios kf^", kf^", fcj'''^, fcg"'" 



PBi < P < PBo 



ji,axis 1 o/l-cdgc r.axis-i 
fi,p -I- (tQ ; 



P ^ PBo 



" Note that the constants ci, C2, C3 and C4 are defined by the values of the free parameters and the conditions that the 
velocity and acceleration are continuous at the two boundaries. 

^ There is a typographical error in the equivalent expression in CL, which should be the same as that given here 
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(a) Total intensity: model 
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(b) Data: counter-jet 


(c) Data: main jet 















Figure 5. A comparison of the model and data in total inten- 
sity at 0.4 arcsec resolution for the irmer ±30 arcsec. (a) Model 
contours. The levels arc logarithmic, increasing by factors of \/2 
from 4.24 - 543 /Li.ly (beam area)~^. (b) Grey-scale of observed 
emission from the coiuiter-jct. The range is — 0.25 mjy (beam 
area) ~ ^ . (c) as in panel (b) , but for the main jet with a rajige of 
0-1 m Jy (beam area) ~ ^ . This display emphasizes the fine-scale 
structure in both jets which cannot be described by the model, 
but shows that the gross features of the predicted and observed 
brightness distributions are very similar. 
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Figure 6. Average transverse profiles of total intensity at a res- 
olution of 2.35 arcsec. The profiles are generated by averaging 
along radii between 45 and 66.5 arcsec from the nucleus and plot- 
ting against angle from the jet axis. The data are represented by 
the full line and the model by the dashed line, (a) main jet; (b) 
counter-jet. The jets are not perfectly straight, so there is a slight 
offset between the observed and model profiles over this range of 
distances. 
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Figure 7. Average transverse profiles of the jet/counter-jet sid- 
edness ratio /j/Zcj at a resolution of 2.35 arcsec. The profiles are 
generated by averaging along radii from the nucleus and plotting 
against angle from the jet axis. The data are represented by the 
full line and the model by the dashed line, (a) 35 — 46 arcsec; (b) 
46.5 - 57.5 arcsec; (c) 58 - 69 arcsec. 
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the latter is nearly independent of distance. [The small offset 
between the observed and model profiles visible in both of 
these figures is caused by slight bends in the jets between 
small and large scales.] 

The following total-intensity features are described ac- 
curately by the model: 

(i) The jets are initially well coUimated and fiare to a 
projected opening angle of ~ 30° at a distance of 20 arcsec 
from the nucleus (Figs|^ and b and^t and b). 

(ii) Both jets are faint and narrow within 5 arcsec of the 
nucleus. 

(iii) Between 5 and 18 arcsec from the nucleus the main 
jet is very bright (Figs|2t ~ c). 

(iv) Further out, the brightness of the main jet decreases 
monotonically, but the longitudinal intensity profile flattens 
with increasing distance (Fig. ^) . 

(v) The counter-jet peaks twice, at 15 and 30 arcsec 
(FigsOK- c). 

(vi) The jet /counter-jet sidedness ratio is high within 
30 arcsec of the nucleus, thereafter maintaining a constant 
value on-axis (Fig. |^. 

(vii) The sidedness ratio is higher on-axis than at the 
edges of the jets (FigsOt, f, h, i; Fig. 0. 

4.2 Linear polarization 

Fig. |H| compares the model and observed degrees of polar- 
ization, p at resolutions of 2.35 arcsec (a - c) and 0.4 arcsec 
(d - f). The values of p are represented as grey-scales and as 
profiles along the jet axis. Fig. [^compares the average trans- 
verse profiles of the degree of polarization at distances from 
the nucleus between 45 and 66.5 arcsec. Finally, Figs llOllllI 
and 1121 represent both the value of p (vector length) and the 
direction of the apparent magnetic field (vector direction) at 
the two resolutions. 

The following polarization features are well described 
by the model: 

(i) The degree of polarization in the brightest portion of 
the main jet has a V-shaped structure, with high polariza- 
tion close to the edges and low polarization on-axis (Figs|HJl 
and e). 

(ii) Both jets show a high degree of polarization at their 
edges, with the apparent magnetic field parallel to the edge 
everywhere (Fig. EH- 

(iii) The main jet shows a transition from longitudinal to 
transverse apparent field on-axis at a distance of «30 arcsec 
from the nucleus whereas the counter-jet shows transverse 
apparent field on-axis everywhere there is adequate signal 
(Fig.Hnj. 

(iv) Within ~3 arcsec of the nucleus, the counter-jet has 
a transverse apparent field with p ~ 0.15 - 0.2. At 3.6 arcsec, 
there is a polarization minimum (consistent with zero); at 
larger distances the degree of polarization increases with 
distance (Fig. |HJ;). The polarization direction is difficult to 
see in Fig. llOf a). where the model is blanked if J < 5a/ 
or P < 3(Tp to match the data. Close to the nucleus, the 
counter-jet has discrete knots of emission with detectable 
polarization; the model is smoother, has a lower peak in- 
tensity and is totally blanked. For this reason Fig. ll2l shows 
model polarization vectors close to the nucleus with minimal 
blanking. 



(v) The counter-jet shows a much more pronounced ridge 
of high, transverse polarization on-axis than does the main 
jet (Fig.iJ. 

(vi) The on-axis polarization of the counter-jet is sys- 
tematically higher than that of the main jet at distances 
>30 arcsec, where both show transverse apparent field 
(Fig.lHl;). 

4.3 Features that are not fitted well 

As discussed in more detail by Worrall et al. (in prepa- 
ration), the bright region of the main jet is resolved into 
complex knots and filaments and its structure is clearly not 
axisymmetric in either total intensity or linear polarization 
(the same is almost certainly true of the counter-jet, but we 
see only its brightest few knots). Our model fits the inner jets 
with a smooth brightness distribution of the correct mean 
value, but there are large local deviations. An inevitable 
consequence is that there are fiuctuations in sidedness ra- 
tio, which we cannot fit. We emphasize that relying on the 
sidedness ratio alone to derive jet properties is dangerous: 
it is necessary to average the intensity over a region large 
enough to contain many fluctuations, but small compared 
with variations in the underlying flow. This may not always 
be possible. The intensity fits close to the nucleus have two 
problems which may be due to small-scale fluctuations: 

(i) There is significant emission from the counter-jet close 
to the nucleus, (Fig.|SJ and the sidedness ratio derived by in- 
tegrating over the emission between 0.8 and 1.5 arcsec from 
the nucleus at 0.4 arcsec resolution is 6.65, much smaller 
than is seen further out (Fig. |3|;) . In principle, such an in- 
crease of sidedness cannot occur in a monotonically deceler- 
ating jet and our model fails to reproduce it. We return to 
this point in Section [5.3.21 

(ii) The first brightening of the counter-jet is slightly, but 
significantly further from the nucleus than the corresponding 
feature in the main jet: the model predicts that they should 
both brighten at the same place (Fig. |3Ji) . This effect causes 
the huge peak in observed sidedness ratio seen in Fig. |2Ig). 
It almost certainly results from the coincidence of a knot in 
the main jet with a minimum in the counter-jet at ~9 arcsec 
from the nucleus (Fig. 01. 

A related problem causes a discrepancy between the pre- 
dicted and observed degrees of polarization observed on-axis 
in the bright region of the main jet: this appears to be due to 
filaments with aligned apparent fields crossing the jet axis at 
an oblique angle (Fig. Illj a'). At high resolution, the model 
fits the lower bound of the polarization profile quite well 
(Fig.H). 

In addition to small-scale, non-axisymmetric structure 
(which affects all of the sources we have studied at some 
level) the model of NGC 315 also fails to fit two other fea- 
tures. 

(i) The on-axis sidedness ratio is slightly higher than pre- 
dicted at large distances from the nucleus, leading to a larger 
contrast in sidedness ratio between centre and edge than is 
predicted by the model (Figs |3 and |7J. This may indicate 
that the transverse velocity profile is more complex than we 
assume, and we discuss this point further in Section f5. 3. 31 

(ii) The average transverse polarization profile for the 
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Figure 8. A comparison of the degree of polarization, p, for model and data, (a) and (b) grey scales of p at 2.35 arcsec resolution in the 
range — 0.7. (c) Profiles of p along the jet axis for the data (full line) and model (dashed). Panels (a) - (c) show the inner 66.5 arcsec 

of each jet. (d) and (c) grcy-scalcs of p for the main jet at 0.4 arcsec resolution, (d) model, (c) data, (f) Profile of p along the axis for the 
main jet. Panels (d) - (f) show the inner 21 arcsec of the main jet. The data are blanked only if / < 5a"/, to avoid discriminating against 
areas of low p. The models are not blanked. 
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Figure 9. Average transverse profiles of the degree of polarization 
for the main and counter-jets at a resolution of 2.35arcsec. The 
profiles are generated by averaging along radii between 45 and 
66.5arcsec from the nucleus and plotting against angle from the 
jet axis. In both panels, the data are represented by full lines and 
the model by dashed lines, (a) main jet, (b) counter-jct. 



main jet is slightly flatter than predicted by the model at dis- 
tances between 45 and 66.5 arcsec from the nucleus (Fig.^J. 



5 PHYSICAL PARAMETERS 

5.1 Summary of parameters 

In this section, all distances in linear units are measured 
in a plane containing the jet axis (i.e. not projected on the 
sky). The parameters of the best-fitting model and their 
approximate uncertainties are given in Table |21 

5.2 Geometry and angle to the line of sight 

The best-fitting model has an angle to the line of sight of 
e = 38° ± 2°, consistent with the range 30° < 6» < 40° esti- 
mated bv lGiov annini ct al. ( 2001) from the core prominence, 
proper- motion measurements and jet/counter-jet intensity 
on pc scales (the last using an isotropic emission model). The 
shape of the outer edge is shown in Fig. I13r a'l . The initially 
well-coUimated jets flare in the inner 15kpc to a maximum 
opening angle «20°, which is maintained over most of the 
rest of the modelled region. The cubic fit to the jet radius 
requires that a conical region should start at a distance of 
49 kpc from the nucleus (90 arcsec in projection on the sky) 
and have an opening angle <6°. This is consistent with the 
observed appearance of the jets at larger distances (Fig.0, 
taking into account the bends at 70 arcsec. NGC315's jets 
are therefore similar to those in other objects in recollimat- 
ing to become almost cylindrical, but only after the slight 
bends which limit our modelling. On scales >500 arcsec in 



projection, there is a second rapid expansion JWilhs et alJ 
19811 lBridlelll982h . 



5.3 Velocity 

5.3.1 Model fits 

The boundaries between the three velocity regimes are plot- 
ted, together with their errors, in Fig. I13f a'). The on-axis 
velocity profile, shown by the full line in Fig. I13f b'l. is fit 
with a constant value of /3 = 0.88 ±0.1 over the inner 7 kpc, 
but this value is determined primarily by the data at dis- 
tances >2.5kpc, where the jet widens and brightens. Be- 
tween 7 and 20 kpc, the jet decelerates uniformly, reaching 
a well-constrained asymptotic speed that is still appreciably 
relativistic {(3 — 0.38 ± 0.03) and persists to the end of the 
modelled region at 38 kpc. The best-fitting velocity at the jet 
edge is slower than the on-axis value, the edge/on-axis ratio 
varying from 0.8 ± 0.2 close to the nucleus to 0.58 ± 0.13 at 
large distances. A constant value of 0.6 - 0.7 would fit quite 
well everywhere, but a constant (top-hat) transverse profile 
is also consistent with the data close to the nucleus. The full 
velocity field is shown as contours in Fig. 1141 



5.3.2 Evidence for acceleration? 

On pc scales, both the jet /counter-jet ratio and the appar- 
ent c omponent speed increase with distance fro m the nu- 
cleus JCotton et alJll999l: iGTovannini et al.ll200j) . although 
no changes in velocity have yet been detected in any indi- 
vidual co mponent. The velocit ies derived from the measure- 
ments in ICotton et al.l |l999) , but taking our best-fitting 
inclination angle of ^ = 37.9° and Ho = 70kms~^ Mpc~^, 
are plotted in Fig. 1151 The velocities from /j//cj are system- 
atically higher than those from component motions at the 
same distance from the nucleus for this Hubble Constant. 
An estimate from Ij/Icj for our 0.4 arcsec image between 0.8 
and 1.5 arcsec from the core in projection (~ 0.6kpc along 
the a xis) is also included. This assumes (as in lCotton et alJ 
Il999^ that the rest-frame emission is isotropic, so that: 



1+13 cos 9 
1 - /3 cos e 



We also plot the results of our model fits for the centre and 
edge of the jet for distances >2.5kpc, where they are well 
determined. Aside from the low value at «0.6kpc, the ve- 
locities are consistent with acceleration from /3 « 0.7 at 
Ipc to /3 ~ 0.95 at 10 pc, an approximately constant veloc- 
ity between 10 pc and 8 kpc and deceleration as described 
earlier. M87 llReid et al.lll98gl: iBiretta. Zhou fc Qwenlll995t 
iBiretta fc Junoil Il995l: lJ unor fc Biretta' '1995*) and Gen A 
iTingav. Preston fc Jaun cev 2001; Hardcastle ct al. 200^ 
also show component speeds which increase from pc to kpc 
scales. Similarly, NGC 6251 has been suggested as an ex- 
ample of an accelerating fl ow because its jet / counter-jet ra- 
tio increases with distance iSudou et al .120001) . although the 
de tection of a counter-je t on pc scales was not confirmed 
bv lJones fc Wehri3 i2002l) . Two possible mechanisms for in- 
crease of velocity on pc scales are thermal accelerati on of 
proton-electron p lasma (iMelia, Liu fc F antuzzd |2002ll and 
magnetic driving JVlahakis fc Konig]||2004l)" 
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(a) Polarization vectors: model 
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(b) Polarization vectors: data 
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Figure 10. Vectors with lengths proportional to the degree of polarization, p, and directions along the apparent magnetic field, su- 
perimposed on grey-scales of total intensity. The resolution is 2.35 arcsec and vectors are plotted every 1.5 arcsec. The polarization and 
angular scales are indicated by the labelled bars in the lower panels and the areas plotted are the same as those in Figs 13 and |H] Vectors 
are plotted only where / > 5a"j and P > Strp (Table 0. (a) model; (b) data. 
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Table 3. Fitted parameters and error estimates. 



Quantity Symbol opt min" max 



Angle to line of sight (degrees) 


e 


37.9 


36.1 


40.7 


Geometry 
















Boundary position (kpc) 


ro 


49, 


,18 


47, 


,37 


50, 


,61 


Jet half-opening angle (degrees) 


?0 


3, 


.75 


0, 


.08 


6, 


,52 


Width of jet at outer boundary (kpc) 


xo 


12, 


.71 


12, 


,10 


13, 


,37 


Velocity 
















Boundary positions (kpc) 
















inner 


Pvi 


7, 


.59 


4, 


.28 


10, 


,10 


outer 




18, 


.07 


16, 


.24 


19, 


,74 


On — axis velocities / c 
















inner 


/3i 


0, 


.88 


0, 


.77 


0, 


,99 


outer 


/3o 


0, 


.38 


0, 


.35 


0, 


,41 


Fractional velocity at edge of jet 
















inner 


VI 


0, 


.79 


0, 


,59 


1, 


,05 


outer 


Vo 


0, 


.58 


0, 


,45 


0, 


,71 


Emissivity 
















Boundary positions (kpc) 
















inner 


Pci 


2, 


.52 


0, 


,00 


3, 


,22 


2 


Po2 


3, 


.53 


2, 


,70 


4, 


,75 


3 


pes 


9, 


.41 


9, 


,14 


9, 


,73 


4 


pBi 


10, 


.05 


9, 


,71 


10, 


,30 


On — axis emissivity exponents 
















inner 


Ein 


3, 


.15 




< 


4, 


,00 


2 


E ■ 


0, 


.03 




< 


2, 


,34 


3 


Emid 


2, 


.79 


2, 


,33 


3, 


,18 


4 




10, 


.39 


6, 


,45 


13, 


,47 


5 


Eout 


2, 


.87 


2, 


,72 


3 


,04 


Fractional emissivity at edge of jet 
















inner boundary 


ei 


1, 


.Ui 


U, 


A 1 




1 o 


boundary 3 


eo 


n 
U, 




n 
u. 




u 




B-ficld 
















Boundary positions (kpc) 
















inner 


PBi 


0, 


.46 


0, 


,00 


6, 


,47 


outer 


PBo 


25, 


.79 


20, 


,85 


31, 


,73 


RMS field ratios 














radial/toroidal 
















inner region axis 


Jl 


1, 


.12 


0, 


,54 


1, 


,71 


inner region edge 


.edge 
H 


0, 


.45 


0, 


,03 


0, 


,82 


outer region axis 


„-axis 


0, 


.61 


0, 


,08 


0, 


,91 


outer region edge 


.edge 
JO 


0, 


.20 


0, 


,00 


0, 


,39 


longitudinal / toroidal 














inner region axis 




1, 


.43 


1, 


,18 


1, 


,70 


inner region edge 


^edge 


0, 


.95 


0, 


,77 


1, 


,14 


outer region axis 


r,axis 


0, 


.97 


0, 


,79 


1, 


,19 


outer region edge 


l.cdge 
Kg 


0, 


.37 


0, 


,19 


0, 


,52 



The Symbol < means that any value smaller than the quoted maximum is allowed. 



As noted bv lCotton et alJ l^99|), it is unclear whether 
the apparent acceleration on pc scales occurs because the 
bulk flow accelerates or because we see different parts of 
a jet stratified in velocity at different distance from the 
core. An acceleration process which only reaches a speed 
of /3 ~ 0.7 at 1 pc is impossible to reconcile with FRI ra- 
dio galaxies being the parent population of TeV blazars, as 
highly relativistic (r > 10 - 20) flow is required on 0.1-pc 



scales in the latter class to produce the observed 7-ray emis- 
sion. Velocity gradients have indeed been inferred on sub-pc 
scales to explain the simultaneous observations of slowly- 
moving rad io components and highly variable TeV e mission 
in blazars (jOhisellini. Tavecchio fc Chiabereel l200.'Th . Simi- 
larly, although the anomalously small velocity estimate at 
0.6 kpc in NGC 315 might well result from a random inten- 
sity fluctuation, 3C 31 (LB) and B2 0326-1-39 (CL) also show 
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(a) Polarization vectors: model 



(b) Polarization vectors: data 



P = 1 



2 arcsec 



Figure 11. Vectors with lengths proportional to the degree of polarization, p, and directions along the apparent magnetic field, super- 
imposed on grey-scales of total intensity. The resolution is 0.4 arcsec FWHM and vectors are plotted every 0.4 arcsec. The polarization 
and angular scales are indicated by the labelled bars in the lower panels and the areas plotted are the same as those in Figs Bl and 151 
Vectors are plotted only where / > Scr/ and P > Scrp (Table 0. (a) model; (b) data. 



lower sidedness ratios close to the nucleus than at the flar- 
ing point, suggesting that the material contributing the bulk 
of the emission in the faint region at the jet base really is 
slower than that further out, where our models are well con- 
strained. We cannot resolve the transverse velocity structure 
of the jets in their innermost regions, so it is possible that 
they have very fast central spines whose emission is Doppler 
dimmed on both sides of the nucleus and that the visible 
emission comes from much slower surface layers. Rapid de- 
celeration of the spine to a speed /3i ~ 0.9 at 2.5 - 3.5 kpc 
from the nucleus would make it visible, dominating the emis- 
sion and causing the sidedness ratio to rise, as suggested for 
3C31 by LB. Further from the nucleus, there are local min- 
ima in the sidedness ratio at 35 and 50 arcsec, separated 
by a maximum at 40 arcsec (Fig. These could be inter- 
preted as changes in velocity in the range 0.32 < /3o < 0.42, 
but they show no coherent pattern, so we believe that they 
probably reflect intrinsic differences between the two jets. 



5. 3. 3 Transverse velocity profile 

The discrepancy between the predicted and observed trans- 
verse profiles of jet/counter-jet sidedness ratio (Fig.|7| pro- 
vides the first evidence that the Gaussian and spine/shear- 
layer velocity profiles we have employed so far may be over- 
simplified. The effect, which is present in all three average 
profiles shown in Fig. |3 is that the sidedness ratio has a 
sharp peak within ~ 5° of the jet axis, drops abruptly at 
~ 10° and has a fiatter wing at the edge of the jet. The 
type of velocity profile which could produce the observed 
results would have a central spine of high velocity {/3 ~ 0.5) 
surrounded in turn by a relatively narrow shear layer and 
an outer wing with /3 ~ 0.2. There are important implica- 
tions for the evolution of the magnetic field, which would be 
sheared only over a narrow range of radii, rather than the 
majority of the jet, and for the physics of the deceleration 
process. The other sources we have studied (LB, CL) show 
no systematic discrepancies, but are less well resolved than 
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Polarization vectors: model 




p = 1 2 arcsec 



Figure 12. Vectors with lengths proportional to the degree of 
polarization, p, and directions along the apparent magnetic field, 
superimposed on grey-scales of total intensity for the model at 
a resolution of 0.4 arcsec FWHM. The polarization and angular 
scales are indicated by the labelled bars and the inner ±8 arcsec 
is shown. Vectors are plotted if / > 2/iJy(beam area)~^ (cf. 
Fig. Illi . This plot shows the model polarization structure of the 
counter-jet near the nucleus, including the null at 3.6 arcsec and 
the region of transverse apparent field closer in. The data are not 
plotted for comparison at this resolution because the signal-to- 
noise ratio in the counter-jet is too low, but these features are 
qualitatively consistent with the low-resolution images (Figs 15b. c 
and ITHb'). 

NGC 315 in regions where they have significant changes in 
sidedness from axis to edge. 



from a smooth extrapolation between smaller and larger dis- 
tances (Fig. 1131 1'). We also resolve the inner boundary of 
this region in 3C 31 and B2 0326-1-39 (LB, CL), modelling 
it as a discontinuous increase in emissivity. The bright re- 
gion comes to an equally abrupt end, marked by an almost 
discontinuous drop in emissivity in NGC 315 (Fig. ll3H ) and 
B2 0326+39 (fig. 16c of CL). B2 1553-^24 may show a sim- 
ilar feature, but is less well resolved (fig. 18c of CL). In all 
three sources, the drop in emissivity occurs well before the 
jets recollimate. 3C 31 does not show any sharp decrease in 
emissivity. 

The phenomenon of sudden brightening and expansion 
at a flari ng point close to the nucleus is very common in FR I 
jets (e.g. iParma et ai~lll987f) . but the equally sudden emis- 
sivity drop at the end of the bright region has only become 
apparent from our modelling of NGC 315 and B2 0326-1-39. 
In both sources, the drop is located just after the start of 
the rapid deceleration (p — p„j ; see Table |3J and coincides 
with it to within the errors. This association reinforces the 
argument that dissipation in FRI jets, leading to enhanced 
radio emission and the production of synchrotron radiation 
at much high er frequencies, is asso ciated mainly with their 
fastest parts jLaing fc Bridlell2004) . 



5.4 Emissivity 

The boundaries between the emissivity regions, with their 
errors, are shown in Fig. I13f c'l and the on-axis profile of 
noB^^" , derived from the emissivity, is plotted as the heavy, 
full line in Fig. ll3f d'). no and B are in SI units. The emissiv- 
ity profile at distances <3.5kpc from the nucleus is poorly 
constrained (Table |5Ji. Again, this is partly because the jets 
are faint and poorly resolved, but more as a consequence of 
the poor fit to the counter-jet (Sections 14.31 1. which causes 
the values to change very little as the model parameters 
are varied. The profile plotted in Fig. I13f d') is essentially de- 
termined by the brightness distribution of the main jet for 
distances <3.5kpc. We model it as an inner region with a 
slope -Bin ~ 3 (<2.5kpc) followed by a transition zone with 
a much flatter slope Briso ~ from 2.5 - 3.5 kpc. The sud- 
den increase in the observed surface brightness of the main 
jet (Fig. 2J is then produced by a rapid expansion of the 
jet at roughly constant velocity and emissivity. The emis- 
sivity profile for the bright, well-resolved sections of the jets 
can be divided into two power-law sections, with indices of 
Smid = 2.8 (3.5 - 9.4 kpc) and Kut = 2.9 (>10.1kpc). 
These are separated by a second short transition zone over 
which the emissivity drops by a factor ~2, modelled as a 
very steep power law with index _Efaii « 10.4 between 9.4 
and 10.1 kpc; this might also be represented as a disconti- 
nuity. The bright region (3.5 - 9.4 kpc) contains complex, 
non-axisymmetric and knotty structure (Figsjljj and lllb 'l. 
whose average is represented well but whose details are not. 
The outer section is, by contrast, relatively smooth. There 
is no evidence for any transverse variation of emissivity at 
the start of the bright region, although this is poorly con- 
strained. In the outer section the edge emissivity is about 
half of its on-axis value (a profile of noB^'^" at the jet edge 
is shown as the light, full line in Fig. ll3t '). 

The bright region is clearly differentiated from the rest 
because its emissivity is a factor of 2 higher than expected 



5.5 Magnetic-field structure 

The variation of magnetic-field structure is illustrated by 
Fig. 1161 which shows grey-scales of the fractional field com- 
ponents: radial, {B^y^'^/B, toroidal, {B'i)^''^ /B and longi- 
tudinal, {BfY^'^ / B as defined in Section 13.3.31 Profiles of 
these quantities as functions of distance along the jet axis, 
z, are plotted in Fig. ll7l for on-axis and edge streamlines. All 
three components vary slowly out to a distance of 25 kpc, 
after which they remain constant. The errors on the radial 
component are much larger than on the other two. On-axis, 
all three components are initially roughly equal (i.e. the field 
is on average isotropic); further out the toroidal component 
increases and the radial and longitudinal components de- 
crease. At the edge of the jet, the toroidal component is al- 
ways the largest, and dominates over the other two at large 
distances (Fig. 1171 1'). The main field components in different 
parts of the jets are shown schematically in Fig. I16ld '). 

It is difficult to assess the uncertainties in the field com- 
ponent ratios from the error estimates in Table |3 so each 
panel of Fig. 1171 includes a shaded area defining the region 
which the profile could occupy if any one of the six free pa- 
rameters defining it is varied up to its quoted error. 

The field structure we infer for the jets in NGC 315 
within «15 kpc of the nucleus shows some sim ilarities to that 
suggested for M87 bv lPerlman et aP 1^9^. In our model 
for NGC 315, the best estimate for the field configuration on- 
axis in this region is roughly isotropic (although there are 
large uncertainties in the magnitude of the radial compo- 
nent; Fig. 117b ): at any rate, the radial component is larger 
on-axis than at the ed ge, where longi t udina l and toroidal 
components dominate. iPerlman et al.l lll999fl suggest that 
the perpendicular component of apparent field in M87 is 
larger close to the axis and is associated preferentially with 
the optical knots whereas the (mainly radio) edge emission 
has its apparent field aligned with the jet axis. In our pic- 
ture, the aligned apparent field at the edge comes from the 




Figure 13. Panels (a) and (c): sketches showing the relative positions of the boundary surfaces between velocity and emissivity regions. 
The boundaries are defined in Tableland their positions for the best fitting model are given in TableEl The full vertical curves indicate the 
boundaries and the shaded areas their allowed ranges, also from Table|3 (a) Velocity. The regions of approximately uniform deceleration 
and asymptotic outer velocity are marked, (c) Emissivity. The region of enhanced emissivity between 3 and 10 kpc is indicated, (b) and 
(d): profiles of intrinsic parameters along the jets in the rest frame, (b) the velocity profile along the jet axis (solid line) and jet edge 
(dashed line), (d) no-B^+" derived from the emissivity, with no and B in SI units. Solid line: model; dashed line, adiabatic approximation 
with the magnetic-field structure expected from flux freezing. The adiabatic model curves are plotted twice, normalized to match the 
model at 3.5 and 30 kpc from the nucleus, respectively. The heavy lines show the on-axis profiles; thin lines the edge profiles. 



projection of intrinsically longitudinal and toroidal compo- 
nents and the radial component is significant only on-axis. 



Finally, we note that symmetry of the average trans- 
verse intensity and polarization profiles rules out a global 
helical field unless the jet is o bserved side -on in the rest 
frame of the emitting material jLainalll98l(l . The condition 
for side-on emission in the rest frame for the approaching 
jet in NGC315, (3 = cos 9 = 0.79, is roughly satisfied before 
the jets decelerate, but not for distances > 10 kpc from the 
nucleus. In the counter-jet, the condition can never hold. 
The intensity and polarization profiles at large distances 
(Figs |S| and |5J are extremely symmetrical, particularly in 
the counter-jet, and we can rule out globally-ordered helical 
fields on these scales. A configuration in which the toroidal 
component is vector-ordered but the longitudinal compo- 
nent has many reversals would give identical brightness and 
polarization distributions to those we calculate for a fully 
disordered, anisotropic field and we cannot rule out such a 
field configuration. 




Figure 14. Contours of the model velocity field. The contours 
are at intervals of 0.05 in the range f3 = 0.25 - 0.85. 
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then shear would be locahsed at intermediate 
radii, a more complex situation than that considered by 
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Figure 15. A comparison of velocity estimates on pc and kpc 
scales, plotted against distance from the nucleus (note the log- 
arithmic scale). We assume that 9 = 37.9° everywhere. Filled 
squares: velocities from proper motions; open squares: velo cities 
from jet/counter-jet ratios (both from ICotton et al.lll99Hl but 
with our choice of Hubble Constant). The open triangle shows the 
velocity derived from the jet/counter-jet ratio at 0.4arcsec reso- 
lution close to the nucleus. All velocity estimates from intensity 
ratios alone are calculated for isotropic emission in the rest frame. 
The full and dotted lines show our model fits for the centre and 
edge of the jets, respectively. They are plotted only for distances 
>2.5kpc, where they are well constrained (see Section l5.3.2l . 



5.6 Flux freezing and adiabatic models 

In this section, we follow lllaing fc Bridlel l^2004^ in referring 
to our detailed fits as free models in order to distinguish 
them from the adiabatic models considered here. Given the 
assumption of flux freezing in a jet without a transverse 
velocity gradient, the magnetic field components evolve ac- 
cording to: 

Br oc {xpry^ 

Bt oc (x/3r)~^ 
Bi (X X 



in the quasi-one-d imensional approxi mation, where x is the 
radius of the jet jBaum et al.lll99iD . The dashed lines in 
Fig. llTl show the predicted evolution of the field components, 
normalized to match the models at a distance of 30 kpc from 
the nucleus. The evolution of the longitudinal and toroidal 
components is qualitatively as expected but quantitatively 
inconsistent: the longitudinal/toroidal ratio decreases with 
distance, but much less rapidly than predicted. Shear will 
tend to slow the decline of the longitudinal component, how- 
ever, so a n axisymmetric adiaba tic model of the type de- 
scribed bv lLaing fc Bridie (|2q3) may provide a better de- 
scription. If the transverse velocity profile indeed has the 
form suggested by the average sidedness profile (Fig. |7| Sec- 



iLaing fc Bridlel 1(20041) . In contrast, the evolution of the on- 
axis radial component, is qualitatively inconsistent with fiux 
freezing in any simple axisymmetric, laminar-flow model: the 
radial/toroidal field ratio decreases with distance instead of 
remaining constant. Anomalous behaviour of the radial com- 
ponent also occurs in 3C31 and B2 03264-29 (LB, CL). 

Assuming that the radiating electrons suffer only adia- 
batic losses, and again adopting the quasi-one-dimensional 
approximation, the emissivity is: 

eoc(a;2/3r)-(i+2"/^)Bi+" 

l|Baum et al.l Il997t iLaing fc Bridlel l2004^ . B can be ex- 
pressed in terms of the parallel-field fraction / = (Bf)^''^ / B 
and the radius x, velocity /3 and L orentz factor F at som e 
starting location using equation 8 of iLaing fc Bridi3 ||2i 



B oc 



1/2 



We compare the adiabatic and free model emissivity profiles, 
normalized at distances of 3.5 and 30 kpc, in Fig. I13f d). 
In the innermost region, where the fit to the counter-jet 
is poor, the slope Ein of the emissivity variation is es- 
sentially unconstrained (Section 15.41 and Table |^. Every- 
where else there is a clear difference: the emissivity pre- 
dicted by the adiabatic approximation falls more rapidly 
than that derived from the free model. Dissipative processes 
must occur to accelerate the radiating particles where X- 
ray synchrotron emission is detected, i.e. at distances up 
to at least ISarcs ec in p rojecti on, correspo nding to 8 kpc 
in the jet fram e jWorrall. Birk inshaw &: H ardcastlj l2003t 
iDonato. Sambrun^^^k)zzill20d4l) . The failure of adiabatic 
models is therefore inevitable at these distances, but our re- 
sults suggest that they are inadequate to describe the emis- 
sivity variations anywhere in the modelled region. The dis- 
crepancies are probably too large to be accounted for by field 
amplification due to shear in a lamin ar, axisymmetric, adia - 
batic model of the type developed bv lLaing fc Bridl3 i2004l) . 
In general, we have found that adiabatic models only come 
close to fitting the brightness distributions of FRI jets after 
they have recoUimated (LB, CL, Laing & Bridle 20041 and 
we have not modelled this region in NGC 315. 



6 SUMMARY AND FURTHER WORK 
6.1 Summary 

We have shown that the synchrotron emission from the flar- 
ing region of the jets in the FRI radio galaxy NGC 315 can 
be fit accurately on the assumption that they are intrin- 
sically symmetrical, axisymmetric, relativistic decelerating 
flows. The functional forms we use to describe the geome- 
try, emissivity, velocity and magnetic-field structure are very 
close to those developed in our previous work (LB, CL). The 
geometry and the relative locations of the emissivity and ve- 
locity regions are very similar to those in two of the other 
sources we have modelled (3C31 and B2 0326-1-39), except 
that all of the physical scales in NGC 315 are larger by a 
factor of «5. 
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Figure 16. Panels (a) - (c): grey-scales of the fractional mag netic-field components, (a) radial, (5^)1/2/^. (b) toroidal {B^y^'^/B; 



(c) longitudinal {Bfy/^/B. B = {B': 
deduced from the model fits. 



■ Bf)^^'^. (d) Sketch of the dominant field components at different locations in the jet, as 



6.1.1 Geometry 

We have modelled only the flaring region within 70 arcsec 
(in projection) of the nucleus, as the jets bend shortly 
thereafter. As in other objects we have studied, the ra- 
dius X of its outer isophote is well fitted by the expression 
X = aiZ^ -\- aaa^, where z is the distance from the nucleus 
along the axis. The jets make an angle oi Q = 38° ± 2° with 
the line of sight, so the size of the region we model is 38 kpc 
and the intrinsic length of the flaring region is «50 kpc, much 
larger than in the other objects. 



6.1.2 Velocity 

The velocity is well constrained from 2.5 kpc outwards, 
where the jet brightens rapidly (the Baring point). From 
2.5 to «8kpc the on-axis speed is consistent with a constant 
value of /3 = 0.88 ±0.11. This is very similar to the values 
derived at the flaring point for the other sources we have 
modelled in detail (LB, CL) and from a statistical analysi s 
of sidedness ratios for a larger sample jLaing et ahlligoi^ . 
An anomalous region of low jet/counter-jet sidedness ratio 
<2.5 kpc from the nucleus appears to indicate a lower veloc- 
ity there, but the jets are faint and poorly resolved, so this 
could be due to local fluctuations in the jet or counter-jet 
brightness. Between 8 and 18 kpc the jet decelerates uni- 
formly to an asymptotic speed of /3 = 0.38 ± 0.03 which is 
maintained until the end of the modelled region. B2 0326-1-39 
and B2 1553-1-24 show velocity profiles of identical form, but 
with lower asymptotic velocities (CL), whereas 3C 31 con- 
tinues to decelerate slowly on larger scales (LB). NGC315 



shows a significant transverse velocity gradient, with an edge 
velocity consistent with 0.6 - 0.7 of the on-axis value every- 
where, as in the other sources. There are hints from the sid- 
edness ratio at large distances that our assumed (Gaussian) 
form for the transverse velocity profile may be inadequate, 
and that a profile with a central spine (/3 ~ 0.5) separated 
from an outer sheath (/3 ~ 0.2) by a relatively narrow shear 
layer may provide a better fit. 



6.1.3 Emissivity 

The emissivity profile along the jets is modelled as three 
main power-law sections with slopes of —3.5 (0 - 2.5 kpc; 
very poorly constrained), —2.8 ± 0.5 (3.5 - 9.4 kpc) and 
—2.9 ±0.2 (9.4 - 38 kpc). These are separated by short tran- 
sition zones, also modelled as power laws. The first of these 
(2.5 - 3.5 kpc) is roughly constant and represents the bright- 
ening of the jet as a very rapid expansion at constant emis- 
sivity. The second, from 9.4 - 10.1 kpc is very steep and 
describes an almost discontinuous drop in emissivity by a 
factor of 2. Simple adiabatic models predict too steep an 
emissivity decline, as we also found for the flaring regions 
of other jets. The emissivity is centre-brightened where its 
transverse variation is well constrained. 



6.1.4 Magnetic field 

To fit the polarization structure of the jets, transverse vari- 
ation of field structure had to be included in our models. On 
axis, the field is roughly isotropic <10 kpc from the nucleus. 
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Figure 17. Profiles of magnetic-field components along the axis of the jet (panels a - c) and its edge (panels d - f). The solid lines show 
the best-fitting model, the shaded areas the uncertainties derived from the limits in Table |^ and the dashed lines the profiles expected 
for a magnetic field frozen into the flow. The profiles for a passively convected field are normalized to the free model predictions at a 
distance of 30 kpc from the nucleus, (a) and (d) toroidal; (b) and (e) radial; (c) and (f ) longitudinal. 



but the radial component declines, leaving an equal mix of 
longitudinal and toroidal field by the end of the modelled re- 
gion. At the edge of the jets, the radial component is small 
and the field configuration evolves from an equal mix of lon- 
gitudinal and toroidal close to the nucleus to almost pure 
toroidal at large distances. From the symmetry of the trans- 
verse intensity and polarization profiles, particularly in the 
outer parts of the jets, we infer that there cannot be a sig- 
nificant, globally-ordered helical field. All three components 
could have many reversals or the toroidal component could 
be globally ordered, provided that the other two are not. The 
evolution of the radial field component along the jets is not 
consistent with fiux freezing in our assumed velocity field. 
That of the toroidal and longitudinal components is quali- 



tatively as expected, but a more complex model, including 
shear, is required for a quantitative test. 

6.2 Further work 

We are currently acquiring VL A data for one further source, 
3C296 ijHarcastle et al.lll997^ . We wiU then present model 
fits for all five sources using the same functional forms in 
order to compare their properties quantitatively. We will 
also investigate more complex transverse velocity profiles, as 
outlined in Section r5.3.3l and develop techniques to deal with 
slightly bent jets. Where the quasi-one-dimensional analysis 
presented here indicates that the adiabatic approximation 
is reasonable, we will fit the bright ness distributions usin g 
the self-consistent adiabatic model o f iLaing fc Bridie J2004 . 
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Once suitable X-ray observations have been made, we plan 
to appl y the conservation-law approach of iLaing fc Bridie 
l^2002b^ to derive the energy and momentum fluxes of the 
modelled jets and their variations of pressure, density and 
entrainment rate with distance from the nucleus. 

There are, as yet, no observations of jets in FRII 
sources, or in any class of source on scales <1 kpc, with res- 
olution and sensitivity adequate for detailed modelling. The 
advent of EVLA, e-MERLIN and broad-band VLBI should 
allow us to apply our techniques to more powerful (and prob- 
ably faster) jets and to probe scales much closer to those on 
which jets are launched. 
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